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Abstract 

Wc study the G/GI/oo queue from two different perspectives in the same heavy-traffic 
regime. First, we represent the dynamics of the system using a measure- valued process that 
keeps track of the age of each customer in the system. Using the continuous-mapping approach 
together with the martingale functional central limit theorem, we obtain fluid and diffusion 
limits for this process in a space of distribution-valued processes. Next, we study a measure- 
valued process that keeps track of the residual service time of each customer in the system. 
In this case, using the functional central limit theorem and the random time change theorem 
together with the continuous-mapping approach, we again obtain fluid and diffusion limits in 
our space of distribution-valued processes. In both cases, we find that our diffusion limits may 
be characterized as distribution-valued Ornstein-Uhlenbeck processes. Further, these diffusion 
limits can be analyzed using standard results from the theory of Markov processes. 



1 Introduction 

Limit theorems for infinite-server queues in heavy-traffic have a rich history starting with the 
seminal paper by Iglehart [20] on the M/M/oo queue. This work initiated a line of research aiming 
to extend Iglehart's results to additional classes of service time distributions. Whitt [-V-i] studies the 
GI /PH/oo queue, having phase- type service-time distributions, and Glynn and Whitt [I I] consider 
the GI /GI /oo queue with service times taking values in a finite set. Furthermore, in [6], [24] and 
[31], the G/GI/oo queue is studied with general service time distributions. [28] gives a survey of 
these results. 

In this paper, we study the G/GI/oo queue as a Markov process. This is accomplished using two 
different methods. In the first method we construct a process that tracks the age of each customer 
in the system and in the second method we construct a process that tracks the residual service 
time of each customer in the system. Although analyzing these processes might at first appear 
to be a complicated task, one of the themes that runs throughout the paper is that techniques 
originally developed for establishing heavy-traffic limits for finite-dimensional state descriptors may 
successfully be applied to the somewhat more abstract infinite-dimensional setting. In our first 
approach we establish fluid and diffusion limits for a measure-valued process tracking the age of 
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each customer in the system using the continuous-mapping approach together with the martingale 
functional central limit theorem. In our second approach we establish fluid and diffusion limits 
for a measure-valued process tracking the residual service time of each customer in the system. 
The representation we use in the second approach was also used by Decreusefond and Moyal [8] 
to analyze the M/G/co queue. Indeed, many of the results and techniques found in this paper 
have been inspired by them. However, our proofs are quite different. In particular, in the second 
approach we establish the fluid and diffusion limits using the functional central limit theorem and 
the random time change theorem together with the continuous mapping approach. We flnd that for 
both the age and residual service time representations the diffusion limit is a distribution-valued 
Ornstein-Uhlenbeck process. We then utilize the highly developed theory of Markov processes in 
order to study our limits. 

Another paper related to ours is Kaspi and Ramanan [2-)]. Although this work analyzes many- 
server queues with general service time distributions, the measure- valued representation of the 
system is similar to our representation. Fluid limits are established for the system in a space of 
Radon-measure-valued processes. However, when establishing diffusion limits for such processes, 
the limit process evidently falls out of the space of Radon-measure-valued process. Indeed, a 
signiflcant challenge in our study was choosing a reasonable inflnite-dimensional space to work in. 
In the work of [8], the space of test functions used is the Schwartz space, or the space of rapidly 
decreasing inflnitely differentiable functions. This space has the disadvantage of not containing test 
functions that would allow one to obtain corresponding heavy-traffic limits for useful functionals 
such as number-in-system and workload. In the present paper, we find that the Sobolev space 
of infinite order (see [KJ], [1]) with respect to L?{^e), where is the excess distribution of the 
service-time distribution, is the tightest space that has all the properties we need to prove limit 
theorems and also enables one to use our results to obtain corresponding limit theorems for useful 
functionals. 

Besides identifying an appropriate infinite-dimensional space to work in, another major contri- 
bution of our work is making a connection between the literature on infinite-dimensional heavy- 
traffic limits for queueing systems ([16], [15], [7], [23], [8], [17], [9]) and the vast literature on 
infinite-dimensional Ornstein-Uhlenbeck processes motivated by applications to interacting parti- 
cle systems ([19], [26], [3], [5], [18], [27], [21], [22], [4]). Our work especially relies on [21] and [22] 
to prove continuity of our regulator map. 

In the forthcoming paper [29], the authors build on the work of [I t] and [21] to prove heavy- 
traffic limits for the G/GI/co queue in a two-parameter function space. They analyze both age 
and residual processes as we do. The main difference between our work and their work is that our 
framework, which uses distribution- valued processes, allows one to apply the continuous mapping 
approach and other standard techniques to obtain the heavy-traffic limits. 

The remainder of the paper is organized as follows. In §2 we derive system equations for both 
the ages of customers in the system and the residual service times of customers in the system. 
These equations will be the starting point for the main results of the paper. In §3 we present a 
regulator map result to be used with the continuous mapping theorem. In §4 we give martingale 
results that will be used with the regulator map of §3, to obtain our fluid and diffusion limits. In 
§5 and §6 we prove our fluid and diffusion limits, respectively. In §7 we analyze the diffusion limit 
for our age process as a Markov process. A corresponding analysis for the limit of residual process 
could be conducted similarly. 
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1.1 Notation 



The set of reals, nonnegative reals and nonpositive reals are denoted by M, and M_, respectively. 
We denote by C°° the set of infinitely difFerentiable functions from M to M and by || • 11^2 the standard 

1 /2 

norm on the space L'^{p), where fi is some measure on M, i.e. for / G L'^i/J-), ||/||l2 = (/r |/P d/i) 
Letting denote the ith derivative of (p E C°° for i > 0, we denote the Sobolev space of order m 
for m > by 

^ 1^ G ||<^|U, ^ IIv^^IIl^ < ooj . (1) 

For each m > 0, this space is known to be a Hilbert space (see §5.2 of [12]). Furthermore, we 
denote the projective limit of the spaces {^^)m>o by 

oo 

$ = Pi q>'n^ 

m=0 

and call <I> a Sobolev space of infinite order (see [10]). It is shown in Lemma 5 of [1] that <I> is 
a nuclear Frechet space with the topology induced by the sequence of seminorms (|| • ||m)m>o- 
Furthermore, since for each m > 0, is a Hilbert space, $ is a countably Hilbertian nuclear space. 
<I> is also a Polish space: It is a complete metric space since it is Frechet and it is separable by 
Assertion 11 of [1]. 

Our primary objects of study are processes that takes values in the topological dual of 
denoted by <I> . To be precise, ^ is the space of all continuous linear functionals on $ and we refer 
to elements of this space as distributions. For G $ and <y9 G $ we denote the duality product of 
/X and if by (/i, 9?) ^ /^(v?)- For G $ , its distributional derivative, denoted by , is the unique 
element of $ such that 

(^/,(^) = 99G^>. 

It is clear that ^ is well-defined by the definition of For G and i G M, we can define TtH as 
the unique element of <I> (when it exists) so that 

where Tt(p is the function defined by Ttip{-) = — t) (when it exists). 

For < T < cxo and Polish space E, we denote by D{\fi,T],E) the space of functions from 
[0,T] to E that are right-continuous with left limits everywhere on (0,r]. We equip this space 
with the Skohorod Ji topology (see [2] or [34]). In the sequel we will be concerned with the cases 
E = M., E = D = D([0, oo),M) and E = ^ . The quadratic covariation of two martingales M 
and N in D([0, oo),M) is denoted by (< M,N >t)t>o ^-nd the quadratic variation of a martingale 
M G L>([0,oo),M) is denoted by (< M >t)t>o = {<~M,M >t)t>o 

In general, nuclear Frechet spaces are infinite dimensional spaces that possess many desir- 
able properties of finite dimensional spaces. For instance, {nt)t>o £ -^([0, oo), $ ) if and only 
if {{iJ,t, f))t>o £ -D([0,oo),M) for each G In proving our main results we make use of the 
following theorem of Mitoma [2o]. 

Theorem 1.1 (Mitoma [2'')]). Let S be a nuclear Frechet space and let (/i")n>i be a sequence of 
elements of D{[0,T],S ), where S denotes the topological dual of S. Then fi"" ^ fx in D{[0,T],S ) 
if and only if for each if ^S, {{fil' , ip))t>o =^ {{nt, ip))t>o in £»([0,r],M). 
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If {fJ.t)t>o £ -C'([0, oo), $') and t G [0, T], we can define the distribution f^fisds as the element 
of ^' such that for ah 99 G 

<||y i^isds,(p^ = j {fis,f) ds. 

Let {J-t)t>o be a filtration on an underlying probability space {il.,J-',¥). A process M € 
-D([0, cxd), <I> ) is a -valued J-t-martingale if for all 99 G ((Mt, (/9))t>o is a M-valued ^j-martingale. 
For two $ -valued martingales, M and in L'([0, cxd), <l> ), their tensor quadratic covariation 
(< M, >t)t>o is given for alH > and all G <^> by 

< M,A^>j (99,1/^) =< (M,V9),(7V.,V) >t, 

and the tensor quadratic variation (<< M >>t)t>o of a ^'-valued martingale M G /^([0, 00), <&') 
is given by (<< M »t)t>o = (< M,M >t)t>o- Two ^ -valued martingales, M and A^, are said 
to be orthogonal if < M,N >= identically. Corresponding notions for the optional quadratic 
variation process [M] are defined analogously. 



2 System Equations 

In this section, we obtain semi-martingale decompositions of the distribution-valued process A = 
{At)t>o, which keeps track of the age of each customer in the system, and the distribution-valued 
process TZ = {TZt)t>o, which keeps track of the residual service time of each customer in the system. 
We begin by treating the age process A in §2.1 and then move on to the residual service time 
process TZ in §2.2. 



2.1 Ages 

Consider a G/GI/ 00 queue with general arrival process {Et)t>o S D([0,oo),M). We denote by 
Ti and rii the arrival time and service time, respectively, of the ith customer to enter the system 
after time 0—, for i > 1. These service times are independent and identically distributed (iid) with 
cumulative distribution function (cdf) F with mean 1, complementary cumulative distribution 
function F = 1 — F, probability density function (pdf) / and hazard rate function h G C^. 

Define {At)t>Q G D{[0,oo), D) so that At{y) denotes the number of customers in the system at 
time t > that have been in the system for less than or equal to y > units of time at time t. At 
time 0—, we assume that there are Ao(y) customers present who have been in the system for less 
than y > units of time and we denote by 

-fi = inf{y > 0\Ao{y) > i} 

the "arrival" time of the ith. initial customer to the system for i > 1. Aq = Aq{oo) denotes the 
total number of customers in the system at time 0—. We also denote by fji the remaining service 
time at time of the ith initial customer in the system. The distribution of fji, conditional on the 
arrival time fj, is given for j; > by 

r ,1 1 — Fi — fi + X) 
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We denote by ff. and hf. the conditional pdf and hazard rate function associated with this distri- 
bution, respectively. 

We now derive system equations for a measure-valued process that tracks the age of each 
customer in service. First note that by first principles we have for y > 0, 

Ao Et 
My) = Yl '^{t-r^<y}'^{t<fj,} + Yl '^{t-r^<y}'^{t-n<v^}■ (3) 

i=l i=l 

Our first result provides an alternative way to write (3). 

Proposition 2.1. For each t > 0, 

Ao Ao{y) Et Et_y 

My) = My) - Yl ^{v^<tA{y+n)} - Y ^{v^>y+r^} + - Y ^ {v.<{t-n) Ay} - Y ^{v.>y}- 

i=l i=Ao{y-t)+l 1=1 i=l 

(4) 

Proof. By (3), we have that 

Ao Et 
My) = Y ^{i-n<y}'^{t<m} + Y '^{t~n<y}'^{t~n<v^} 

i=l i=l 

Ao{y) Et 

= Y '^{t~n<y}'^{t<v,} + Y■^{^~^^<y}'^{'^~^^<v^} 
i=l i=l 

Aoiv) Et 
= My) + Y (ht~n<y}ht<fi,} - 1) + + Y^Mt-n<y}ht-n<v.} - (5) 

i=l i=l 

However, 

^{t-n<y}'^{t<ri,} - 1 = ^{t-n<y}'^{fi,<t} + ^{t~n>y} (6) 
= {'^{t-f^<y}'^{m<t} + '^{t-f^>y}'^{-f^+iii<y}) + ^{t-n>y}'^{-f,+i),>y}- 

and, similarly, 

^{t-n<y}'^{t-n<r,,} " 1 = '^{t~n<y}'^{rj,<t~T,} + '^{t~n>y} (7) 

= {'^{t-r,<y}'^{'m<t-Ti} + l{t-r,>s/}l{»?i<j/}) + '^{t-T,>y}'^{r],>y} , 

Substituting (7) and (6) into (5) and summing over Ao(y) and Et completes the proof. □ 

We now provide an intuitive explanation for the terms appearing in (4) . The first term represents 
the number of customers in the system at time 0— that have been in the system for less than or 
equal to y units of time; the second term represents the number of departures by time t of those 
initial customers that had total service times less than or equal to y units of time; and the third 
term represents the number of initial customers that had been in the system for less than or equal 
to y units of time at time 0— but have been in the system for time greater than y units of time at 
time t. The fourth, fifth and sixth terms represent similar quantities but for customers that arrive 
to the system after time 0—. 
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Next, for t,y > 0, define = {D^)t>o G D{[0,oo),D) by 

^0 / r^iAtA{y+fi) 



^ / rVi/\t/\(y+Ti) \ 



i=l 

and define D = (A)t>o e -^([0, oo),D) by 

-Et / friiA{t-Ti)Ay 



A(?/) = 2^ y-{Vi<{t-n)Ay} - h{u) duj . (9) 

It then follows from (4) that 

rVi^tHv+n) ^ (■rh'\{t-Ti)Ay 

At{y) = ^o(y) + Et- - Dt{y) - V / hf^{u) du-Y^ h{u) du 

My) Et-y 

'^{-h+m>y} - l{r?.>s/}- (10) 

i=Ao{y-t)+l i=l 

To study the age process as a distribution-valued process, we use a Sobolev space of infinite 
order as the test function space, where $_4 is defined as in (1) with ^ set to the excess of the 
service time distribution: 

^l{A) = [ F{y) dy, for A G B(R+) (11) 
J A 

Since h £ C^, for each tp G the integrals {J^, (p) and {J^e, 9?) both exist, where JTg denotes the 
distribution associated with F^, the cdf of the stationary excess distribution of F, i.e. Fe{y) = 
J^F{x)dx. 

We associate with the process A defined above the process A taking values in <I>^, as defined 
in §1.1, such that for each t > and 93 G <l>_4, 

{Au^) = [ ip{y)dAt{y). (12) 

JR 

We similarly define corresponding processes, and V, associated with and D, respectively. 
We also associate with ^0 a ^^^-valued random variable ^o- It is easy to see that for each t > 0, the 
quantities At, and Vt are well defined elements of It is also easy to see by right-continuity 
of the sample paths of A, and D that G L>([0, 00), ^>^). 

By integrating test functions ip G with respect to each of the terms in (10) it follows that 

^0 r-n+imr^t) ^* m/^it-n) 



ipiy)h{y)dy-y2 ^{y)h{y)dy 
i=i ^0 

. / My) Et-y \ 

- / ^(y) d Yl ^{-n+m>y} + Y ^{^^>y} ■ 

■^^+ \i=Ao{y-t)+l i=l / 

We then have the following two propositions, which allow us to simplify (13): 
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Proposition 2.2. For each t > 0, 

^0 /■- T,+(^,At) [■m'\{t~n) 



^{y)h{y) dy + yZ fiy)Ky) dy = {As, ^h) ds. 

-T, ,_i JO Jo 



i=l 1 = 1 

Proof. 

^0 r-n+iViM) ^n^A[t-n) 



X] / ^{y)h{y)dy + Y1 / v>{y)Ky)dy 

J—Ti JO 

Y] / 1{0<^<7?JV('5 - n)h{s -fi)ds + y] / l|o<,_^^<^.}V'(s - Ti)h{s - n) ds 

/ 5^ 1{0<S<^JV'(S - fi)h{s - Tj) + V l{o<s-r,<r,,}¥'(s " Ti)h{s - Ti) ds 
■^0 V*=l ^=l J 



{As, fh) ds. 



Proposition 2.3. For each t > 0, 

Ao{y) Et 



□ 



/ ^{y)d\ '^{~n+m>y}+^^H>y}] =EMO)+ [ {As,(p')ds. 



j{y~t) 

Proof. Integrating by parts, we have that 

A)(y) Et-y 

^{vi>y} 



•^^+ \i=Ao{y-t)+l i=l 

r. I Ao(y) Et-y 

= / Y M-h+n,>y} + Y '^{m>y} V^'l?/) dy 

''^+ \i=Ao{y-t)+l i=l ) 

= / ^{r^>-y-T,+f,i>y,h+y<t} + Y ^{rh>y,n+y<t} f'{y) dy 

'^^+ \i=l i=l J 

Ao „ Et „ 

= S/ '^{n>-y-h+fii>y,n+y<t}^'{y)dy + Y '^{m>y,ri+y<t}V\y) dy 

^0 r-t Et ,.t 

i=i i=i 

.t / Ao ^ Et \ 

= / (Yl ^{0<s~n<~n+ij,}V>' (S -Ti) + Y^ l{0<.-r,<r,J'/''(s - n) ds 
■^^ \i=l 1=1 / 



ip {u)Qs{du)j ds 
7 



(^As,ip'^ ds. 

□ 

Combining Propositions 2.2 and 2.3 with system equation (13), we arrive at 

{At, V?) = {Ao, V?) + {£t - - Vt, if) - f {As, hip) ds+ [ {As, if') ds, (14) 

Jo Jo 

where we define the $^-valued process £ = E.5q so that {£t,^) = Etip{0) for each ip G <I>^ and 
t > 0. In general, we refer to (14) as the semi-martingale decomposition of A for reasons that 
become clear in §4. 




2.2 Residuals 

We next move on to the residual service time process TZ. As in §2, customers arrive to the system 
according to a general arrival process {Et)t>o £ -D([0, oo),M) and we denote by Tj and 77^ the arrival 
time and service time, respectively, of the ith customer to arrive to the system after time 0— . 
Customer service times are iid with cdf F. We assume the service-time distribution has a bounded 
hazard rate function, but here we make no assumptions on the smoothness of the hazard rate 
function as we did in §2.1. Assuming the boundedness of the hazard rate is helpful in defining our 
space of test functions. In general, we could drop all assumptions on the service time distribution, 
but this would require us to restrict our space of test functions. 

Let Rt{y) denote the number of customers at time t > that have less than or equal to y G M 
units of service time remaining. Note that as in [8] , we allow y < so that in addition to keeping 
track of customers present in the system at time t, we also keep track of customers who have already 
departed the system. We assume that at time 0— there are Ro{y) customers in the system that 
have less than or equal to y > units of service time remaining. We let Rq = i?o(oo) denote the 
total number of customers present in the system at time 0—. By first principles, it then follows 
that 

Et 

Rt{y) = Ro{t + y)+Y. Hri,-{t-n)<y}- (15) 

i=l 

The following Proposition presents an alternative form of (15). 
Proposition 2.4. For each t > and y G M, 

Et Et 

Rt{y) = Roiy) + {Ro{t + y)- Ro{y)) + ^ '^{m<y} + ^ '^{m>y,{n+n.)-t<y}- (16) 

i=l i=l 

Proof. By (15), 

Et Et 

Rt{y) = Ro{t + y) + Yl Mv.-{t-n)<y} = ^0(2/) + {Ro{t + y)- Roiy)) + Yl hm-it-n)<y}- (17) 

i=l i=l 

However, 

l{r,,-{t-r,)<j/} = l{r,,<y} + '^{v.>y,V^~{t~n)<y}■ (18) 

Substituting (18) into (17) and summing over Et, completes the proof. □ 
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We now give an explanation for each of the terms appearing in (16). The first term represents 
the number of customers in the system at time 0— with less than or equal to y units of service time 
remaining. The second term represents the number of customers that arrived to the system with 
greater than y units of total service time but at time t have less than or equal to y units of service 
time remaining. The third and fourth terms can be explained analogously but for customers that 
arrive to the system after time 0—. 

For each y > and t > 0, define G = {Gt)t>o G D{[0, oo),D) by 

Et 

Gt{y) = ^ (1|,,<,} - F{y)) . 

i=l 

By Proposition (2.4), Rt{y) may then be rewritten as 

Et 

Rt{y) = Ro{y) + (i?o(i + y)- Roiv)) + Gtiy) + EtF{y) + ^ ^n.>y,{n+v.)-t<y}- (19) 

1=1 

To study this residual process as a distribution-valued process, we will use another Sobolev 
space of infinite order $7^ as the test function space, where $7^ here is defined as in (1) with 
/i = ^_ + where 

M+(A) = / F(y)dy, for AeB{R+), (20) 
J A 

and fX- is Lebesgue measure on ]R_. Notice /i+ here is defined just as /i in (11). The assumption 
that h is bounded implies that for each ip £ $7^ the integrals {J-, if) and {J-e, ^p) both exist. 

Now associating $^-valued processes TZ, Q and J- to R, G and -F, respectively, as in (12), and 
plugging in test functions and integrating each of the terms in (19), we get that for each (p G ^>7^, 

(7^t, if) = (7^o, ip)+ f ^(y) diRoit + y)- Ro{y)) 

+ {Gt, if) + Et{T, (f) + ip{y) d V ^{^,>y,(ri+r,,)-t<y} ■ (21) 

Jr / 

The following proposition now allows us to simplify the form of (21): 
Proposition 2.5. For each t > 0, 

[ ifiy) d{Ro{t + y)- Ro{y)) + / ip{y) d (j^ Hv.>y,in+m)-t<y}] = " / i^^s, ds. (22) 

JR JR \i=l J "'0 

Proof. The proof parallels the proof of Proposition 2.3. Integrating by parts, we have that 

Et 



i)-t<y} 



/ ^p{y) d{Ro{t + y) - Ro{y)) - / ^{v) d \y2^H>y,{n+vi 
JR JR \i=l_ 

^'^{y<m<t+y} \ f'iy)dy + ( l{r,,>j/,(r. 



+v^)~t<y} ^ [y) dy 
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Ro „ Et „ 

Yl / ^{y<v^<t+y}'P'(y)'^y + Yl / ^H>y,{n+v.)-t<y}V''iy)dy 

i=l r=l 
Ro f Et „ 

/ '^{m-t<y<mW iy) + Y / '^{{r.+m)-t<y<mW iy) "^y 

Ro rt Et ,.t 



y] f ^'(fii- ^)ds+y] f '^{r,<s}^' in + vi - S) ds 

i=l -^0 i=l -^0 

i f Rs{du) 1 ds 
Jo \Jr+ J 



{TZs,<p ) ds. 







□ 



Substituting (22) into (21), 

{nt,ip) = {no,ip) + {gt,^) +Et{j^,^) - [ {ns,^)ds. (23) 

Jo 

We refer to (23) as the semi-martingale decomposition of TZ. In §4, we will prove that the process 
Q in (23) is a martingale. Note the similarity of (23) with (4) of [8]. 

3 Regulator Map Result 

In this section we show that given a Sobolev space of infinite order $, the integral equation asso- 
ciated with (i^o,Ai) e X £'([0,00),$'), 

{iyt,(p) = {uo,(p) + ifj-t,^) + {vs-,Bip)ds, (24) 

Jo 

for B a continuous linear operator on $ and if £ defines a continuous function : $ x 
Z)([0, 00),$ ) Z)([0, 00),$ ) mapping (i^OiA*) to z/, under some mild restrictions on B. 
Before we prove our result we need the following definition from [21]: 

Definition 3.1. A family {St)t>o of linear operators on <I> is said to be a (Co, 1) semi-group if 

1. So = I, where I is the identity operator, and for all s,t>0, SgSt = Ss+t- 

2. The map t^Stf is continuous for each ip £ ^. 

3. For each q >0 there exist numbers Mg, cTq and p > q such that 

\\SM\g < Mge'"'^\\^\\p, 

for all (f £ ^, t> 0. 
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Theorem 3.2. Let B he the infinitesimal generator of a (Co, 1) semi-group {St)t>o- Then for each 
(fo,/^) G $ X D([0,oo),<^ ), the equation (24) has a unique solution given by 

{ut,ip) = {i/o,St(p) + {nt,v) + {fJ.s,St^sBip) ds. (25) 

Jo 

Furthermore, (24) defines a continuous function '■ ^' x -^([0, oo), $') — > D([0, oo), $') mapping 

(l/Q,/^) to V. 

Proof. That is a well-defined function from -D([0, cxd), <!>') to -D([0, oo), $') and the form of the 
solution (25) follow from Steps 1-3 of the proof of Theorem 2.1 of [21] (see also Corollary 2.2). 

To show continuity we adapt the argument in the proof of Proposition 3 of [n] (see also 
[22]). By the form of (25), it suffices to show that for each T > the function mapping $ to 
Z)([0, oo), $ ) defined by z/q •— > S*vq and the function mapping L»([0,oo),$ ) to L>([0,r],R) defined 
by /i /o B*S*_gfis ds, where B* and denote the adjoint operators of B and 5^, respectively, 
are continuous. The desired result then follows from Theorem 1.1 and the fact that the addition 
map on D([0, r],]R) x L'([0, r],]R) is continuous at continuous limits (Theorem 4.1 of [32]). 

Let (fo)n>i be a sequence in $ converging to vq. Then by Proposition 0.6.7 of [;]0], for each 
precompact set K C ^ we have 

sup |(z/q — v)l ^ &s n ^ oo. (26) 

Since the map t — > St(p is continuous for each (p £ ^ hy the definition of (Co, 1) semigroup and 
[0,T] is compact in R, the set {SuBip, u G [0,T]} is compact in Thus, applying (26) to the 
compact set K = {5'„i?(/9, u G [0,T]} gives us 

sup \{vq — 1^0, Su^)\ — > as n — > oo. 
ne[o,r] 

By Theorem 1.1, this proves S*u^ =^ S*v in D([0,T],]R) as n ^ oo and thus that the map 
z^o ^ S*vq is continuous. 

Now let {fi^)n>i be a sequence in D([0,oo),<I> ) converging to fi. Then there exist increasing 
homeomorphisms (A")n>i of the interval [0, T] such that for each if £ ^ 

||(^//" — ^;^n(.), (/j^lly — > and II A" — ello" ^ as n — > oo, 

where e denotes the identity mapping on [0, T]. Again, by Proposition 0.6.7 of [30], for each pre- 
compact set K C ^, we have 

sup sup - ^A"(t),¥')| ^ as n^oo, (27) 

tG[0,T] ipeK 

Just as above, the set {SuBip, u £ [0, T]} is compact in Thus, applying (27) to the compact set 
K = {SuBip, u £ [0, T]} gives us 

sup sup — /U;^n((), 5ui?(/3)| — > as n — > oo. 

ie[0,T] ne[0,T] 
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Therefore, 



sup 

t6[o,r] 



sup 

te[o,T] 



< sup 

tG[0,T] ^0 



/ (/^" - tJ'\"{s),St-sBip) ds 
Jo 

Jo 



< sup / Kij,"^ - nx"{s),SuBip)\ ds 

t,u£[0,T] Jo 

<T sup Kn"^ - fixn(^s),SuBif)\ ^ 0, 

s,u€[0,T] 

as n — > oo. By Theorem 1.1, this proves that the map fj. i— > Jq B* S*_gi-is ds is continuous and 
completes our proof. □ 

3.1 Ages 

If we define the hnear operator B-^ on <&_4 so that B-^(p = ip' — hip for ip G $_4, we can write (14) as 

^ = ^'B^(A,'f -^°-^)- (28) 

We now verify that B"^ generates a (Co, 1) semi-group so that Theorem 3.2 will apply to (14). 

Proposition 3.3. B"^ generates a (Co, 1) semi-group {Sf')t>o defined by 

Sf^ = F-\^t{F^) for ^e^A- (29) 

Proof. First we check that B-^ is the infinitesimal generator of the semigroup given by (29). For 
all (/3 G <I>^ we have 



lim 

h^O h 



li^ F-^F{. + h)^{. + h)-^ ^ F{. + h) ^{. + h)-F^ 

h~>0 h 

F-1 [Fip] = F-^ (Pip' - fip) =ip' -hip = B^tp 



h 



Now we check that the semigroup {Sf-)t>o is a (Co, 1) semigroup. Part 1. of the definition of 
(Co, 1) semigroup is clearly satisfied. Part 2. follows from Lemma 2 of [8]. For Part 3., by the 
definition of the seminorms inducing the topology on (1), it is enough to show that for each 
n > there exists an M„ such that for each ip € and s > 0, 



(30) 



i=0 
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Using the chain rule and triangle inequality, for each n > we have 



E 

4 = 



i=0 
n 

= E 

1=0 
n 

i=0 



F 



n+t) 



(n-i) 



F 

{h-h{- + t)Y 



L2 



L2 



L2 



2\\h\\ 



F{. + t) 



F 



L2 



where || • ||oo denotes the sup norm. Now focusing on the norm in the expression above, 



F{- + t) 



F 



L2 



F{x + t) 
F{x) 



+ ) F{x)dx 



lk±Jl^(i){x + tfF{x + t)dx 
F{x) 

< I ip^'\x + tfF{x + t)dx 



so that finally we have (30) with M„ = maxo<i<„ { (")2||/i||J^-*}. 



□ 



3.2 Residuals 

If we define the linear operator on $7^ so that B'^ip = —93 for 93 S $7^, then we can write (23) 
as 

n = ^'^(no,g + Ej^), (31) 

We now verify that B^ generates a (Co, 1) semi-group so that Theorem 3.2 applies to (23). 
Proposition 3.4. B^ generates the (Co, 1) semi-group {Tt)t>o- 

Proof. First we check that B^ generates the semigroup {Tt)t>o- For each ip E we have 



Thip-(p ip{--h)-ip 
lim : = lim - 



h^O h 



h 



-ip 



We now check that {Tt)t>o is (Co, 1) semigroup. It clearly satisfies Part 1. of the definition of 
(Co, 1) semigroup. Part 2. again follows from Lemma 2 of [n]. To show Part 3. note that it suffices 
to show that for all n > 1 
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Then, we have 



|^H(. -t)||22 = / <^i^)(x-tfdx+ ^^'^\x-tfF{x)dx. 

< [ if^''\xfdx+ [ ip'^''\xfF{x + t)dx. 



< / / ^'^''\xfF{x)dx. 



L2- 



□ 



4 Martingale Results 

In this section we show that the processes + P and g defined in §2 are <l>^ and <I>^-valued 
martingales, respectively. That fact that + P is a martingale will be used with the martingale 
functional central limit theorem and the continuous mapping theorem in §5.1 and §6.1 to prove 
fluid and diffusion limits, respectively, for our age process. The fact that ^ is a martingale will 
not be needed to prove limit theorems for our residuals process but can be used to show that the 
diffusion limit for our residuals process is a Markov process (see §7) and possibly in other future 
work. 

4.1 Ages 

First we first show that the process V + defined in §2.1 is a martingale with respect to the 
filtration {J-'j^)t>o defined by 

.Fj-^ = a{l{^^=o}, l{r,^<s-n},s < t, i = 1, 2, Et} V a{Es, s < t} V AA. 

Proposition 4.1. The process V^+T) is a ^ j^-valued Tf" -martingale with tensor quadratic variation 
process given for all ifjip £ by 

+ V »t{Lp,il)) I ip{x-fi)'tp{x-fi)hf^{x)dx + 'S2 ip{x)^p{x)h{x)dx. 

i=i i=i 

(32) 

Proof. We first analyze V. Note that by (9) for t > 0, we have that 

oo 

where for each i > 1 

rrjih{t-Ti)hy 

Dliy) = i{r,,<{t~n)+Ay} - h{u)du. (33) 

Jo 

We will show for the associated valued processes 2?*, z > 1: 
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1. For each i > 1, is a ^>^-valued ^^-^-martingale. 

2. For i ^ j, D"^ and are orthogonal. 

3. For each i > 1, the tensor quadratic variation of is given for all ip, ip G <1>_4 by 

10 



ip{x)'4}{x)h{x) dx. 



If we can then show that "^^^=1 {^^h ^) dominated by an integrable random variable uniformly 
over > 0, it will then follow by Lesbegue's dominated convergence theorem for conditional 
expectations and 1. above that 



E [{Vt, if) \J^s] = E 



i=l 



i=l 



1=1 



and hence T> will be a valued ^j-martingale. Note that for each /c > 1, we have that 



1=1 



E 

1=1 

oo 



(p{x)d[ l{rn<x} 



h{u) du 



sup \if{s)\ 1{ <(t_^.)+|+ / h{u)du\ 

j=l o<s<t \ Jo J 

<Et sup |(^(s)|(l + t||/l||oo). 

0<s<i 

However, clearly K[EtsnpQ^s<t\ip{s)\{l + t\\h\\oo)\J^s] = supo<s<t + t\\h\\oo)^Et\J^s] < oo, 
which completes the proof that 2? is a <l>^-valued .F^-martingale. 

Further, a similar dominated convergence argument along with 2. and 3. above shows that the 

quadratic variation of the $_^-valued ^^-^-martingale 2? is given for all 93, £ by 

« V. »t (99, tp) = '^« V' »t {if, tp)=^ I ip{x)ilj{x)h{x) dx. (34) 



i=l 



i=l 



It remains to show 1.-3. We begin with 1. It suffices to show that for each (p E S, ((fj, (p))t>o is a 
real-valued J^^^- martingale. First, {{Vl, (p))t>o is clearly .F/^-adapted. To show that the martingale 
property holds for {{V^, V3))t>o, it suffices to show that the martingale property holds for D^{y) for 
each fixed y > 0. It will then follow that for s < t, 



Dl{y)^{y)dy\Tf 
E[Dl{y)\J^f] ^{y)dy 
Di{y)^{y)dy, 
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and so will be a <I>^-valued J^^-martingale. Since y f\{t — t,i) = {t f\ (xj + y) — Tj), we have 

We know by the proof of Lemma 3.5 of [24] that D*(oo) is an .Fj^-martingale and it is easy to 
see that + y is a .7^/^-stopping time. Therefore, the stopped process (D*^^^ _^^j(oo))t>o, is a 
.F/'-martingale. 

We now focus on 2. To prove orthogonality of and T)^ for i 7^ j, it suffices to show that 
for all ifjip & the process {{Vl,ip){Vl,'ip))t>o is a .F/^-martingale. As in 1., it suffices to prove 
that for each fixed y > the martingale property holds for the process {Dl{y)Dj{y))t>o. Again, 
this follows from the fact that {Dl{oo)D'l{oo))t>o is a martingale, which we know from the proof 
of Lemma 3.5 of [24]. 

We now calculate the tensor-quadratic variation for to prove 3. First by (33), for all (p G 
we have 

{Vl,^)= ^{x)dDi{x)= ' 99(x) d ( 1|^^<,} - / ' h{u)du]. (35) 

JM+ Jo \ Jo / 

Therefore, as on page 259 of [24], 

«{Vl^p)»t=J if{x)d{l{rn<x}- j h{u)duj=j if{xfh{x)dx, (36) 

so that 

= - (« (p?, V? + ^) »t - « {v\ V9 - V) »t) 



1 / r^i/\{t~n)+ rv^A{t-Ti)+ 



4 

/ /■»?iA(i-T,)+ ^ 

{^p{x) + 'il){x)Y h{x) dx — I {ip{x) — 'ilj{x))'^h{x) dx 



4 



ip{x)'}jj{x)h{x) dx 

where the second equality follows from polarization and the third equality follows from (36). 
We now analyze . First note that 

A°(y) = EA°"(y)' (37) 

i=l 



where for each i > 1, 



Dt'\y) = ^{ij,<tA(y+n)} - / hf^{u)du. 



Replicating the analysis for "D above, we can show that for each i > 1, T>^'^ = ('D^'*)j>o is a 
valued martingale and that for i ^ j, T)^'^ and V^'^ are orthogonal. We now calculate the 
tensor-quadratic variation for D^. First, for all (p € we have 

(2?°^ [ ^{x) dD'^/ix) = f ifix - n) d f 1|^^<,| - r''^ hf^{u)du) . (38) 

Jr+ Jo \ Jo J 
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This gives us 

i f^'^^ 

«{T>/,(f)»= / (p'^{x - fi)hfi{x)dx, 
Jo 

Using the polarization identity as in the analysis of P, it then follows that 

« D^'' »t {ip,ip) = ip{x- fi)^p{x - fi)hf^{x)dx. 

Jo 

Summing the quadratic variations of each of the terms in (37) and noting the orthogonality of 
the martingales in the sum, 

((/j,V)= V«2?o.*»t(y^,V') = y" / ipix-n)i;{x-n)hf^{x)dx. (39) 
i=i ^=l -^0 

Now since V and P° are both ^^^-martingales, P'^ + P is a ^^'^-martingale. Furthermore, P 
and are orthogonal since they are independent. Therefore, << + P >>=<< >> + << 
V ». Plugging (34) and (39) into this equality gives us (32). □ 



4.2 Residuals 

We now show that the process Q defined in §2.2 is a martingale. This will be useful in future work 
where we wish to show that the residual service time process is a martingale. Let J-f be the natural 
filtration generated by Q. We then have the following result. 

Proposition 4.2. The process Q is a ^-j^-valued -martingale with tensor optional quadratic 
variation process given for all ^fjip & by 

Et 

[g]ti^,^l;) = ^ipirJ^mV^). (40) 

i=l 

Proof. We first prove the martingale property. Define the filtration {T~ik)k>i by 7ik = a{Et,t > 
0} V a{rii,r]2, .-.,%} V AA. Define the discrete-time D-valued process (G*)j>i by 

k 

G'=(y) = ^ (!{,,<,} -F(y)), (41) 

i=l 

and let {G'')k>i be the associated $7e-valued process. It is then clear by the independence of the 
service times from the arrival process that for each if G $7^, the process {{G'' ,^))k>i is an Tlk- 
martingale. However, since for each t > we have that Et is a stopping time with respect to 
the filtration {'Hk)k>i, it follows that the filtration iJ~iEt)t>o is well-defined and, furthermore, by 
the optional sampling theorem, we have that for each ip S ^ti, {{Gt,'p))t>o = {{G^\<^))t>o is a 
7YEt-™artingale. The result now follows since any martingale is a martingale relative to its natural 
filtration. 

The form of the tensor optional quadratic variation (40) is immediate by Theorem 3.3 of [is]. □ 
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5 Fluid Limits 



In this section, we begin proving our weak convergence results. We consider a sequence of G/GI /oo 
queues indexed by n > 1, each following the assumptions of §2. We assume that the service time 
distribution is held fixed across the systems. We add a superscript n > 1 to all processes and 
quantities defined in §2 to indicate association to the nth queue in the sequence. We focus on fluid 
limits for the age and residual processes in §5.1 and §5.2, respectively. We move on to diffusion 
limits in §6. 

5.1 Ages 

We start with the age processes of §2.1. Define 

An _ _ _ T)" _ T)0,n 

A'' = —, = ^" = — , P" = — , pO'" = i_, (42) 

n n n n n 

and = E'^6o for n > 1. Then by (28) for n > 1 we have 

^" = ^^^(^[J,^" -pO'"-P"). 

We now prove convergence of V^'^ + jointly with {Aq, E^): 
Proposition 5.1. // {Aq,E'^) =^ {Ao,E) in <I>^ x D as n ^ oo, then 

(^[J,^",P°'" + P") ^ (A,'?,0) in X D([0,oo),«>')2 as n ^ oo. 
Proof. We first prove 

pO.n^pn^Q .^^ D([0,oo),^>') as u ^ OO. (43) 
For each (p,ip ^ <I>_4, T > and < t < T we have 



V / )^(x - ff )/ifn (x) dx + \ ip{x)il^{x)h{x) dx 

< — ^ |> / \(p(x)'ip(x)h(x)\ dx +} / \ip(x)ilj(x)h(x)\ dx 



^=0"" i=0 

<-{A^ + E^)T sup \^is)\ sup 1^^(3)111/11100^0, 

n 0<s<T 0<s<T 

in D as n ^ oo by the assumed convergence of Aq and -E". Therefore, (43) follows by the martingale 
FCLT. Then joint convergence holds by virtue of Theorem 11.4.5 of [34] and the fact that the limit 
in (43) is deterministic. □ 

We then have 
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Theorem 5.2. {Ao,E) in^^x D as n ^ oo, then 

^ A in D([0,oo),$') as n^oo, 
where A satisfies the deterministic integral equation 

{At,^) = {Ao,^) + EtifiO) - [ {As,hip)ds+ [ {As,v')ds, (44) 

Jo Jo 

for all if G <I>^. 

Proof. By the assumption and Proposition 5.1, we have 

+ ^ (A,'f,0) in X Z)([0,oo),«>')2 as n^oo. 

Then, since A"" = ^^^(^q,i5" — P'''" — P") and ^'^^t is continuous by Theorem 3.2 and Proposition 
3.3, the resuh fohows from the continuous mapping theorem (see [2] and [34]) with ^'_4 and the 
addition map. The addition map is convergence preserving here by Theorem 4.1 of [32] since the 
Umits of T)^''^ and P" as n — > oo are continuous. □ 

Remark 5.3. Note that we can now use Theorem 5.2 along with Theorem 3.2 to write an explicit 
expression for A. Similarly, we can get explicit expressions for TZ, A and TZ in Theorems 5.5, 6.5 
and 6.7 below. 

5.2 Residuals 

We now proceed to prove a fluid limit for the residual processes of §2.2. Let 

7?" 7?" 

n n n 

for n > 1. Then by (23) we have 

^^ = ^-^^(^«,g"+^"^). 

We first prove convergence of Q"' jointly with (TZq,E"'): 
Proposition 5.4. If{'RJ^,E'') =^ {Ho.E) in^^xD as n ^ oo, then 

(:^[},<f",g") ^ (:^o,'?,0) in X D([0,oo),$')2 as n ^ oo. (45) 
Proof. Notice that for each if G {Q^, ip) can be written as 




The first term above converges to the function by the functional weak law of large numbers. Thus, 
by continuity of the composition map at continuous limit points (see page 145 of [2] or Theorem 
13.2.1 of [34]), we have {Q^,ip) in as n ^ oo. Since this limit is deterministic, we get the 
full joint convergence (45) by Theorem 11.4.5 of [34]. □ 
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We then have 

Theorem 5.5. // {H}^, {Hq, E) in <^'^ x D as n oo, then 

'RJ'^IZ in L»([0,oo),^'') as n ^ oo, 
where TZ satisfies the deterministic integral equation 



{Tit, = {Tlo, + Et{T, if) - [ {n,,^) ds, (46) 

Jo 

for all if £ ^Ti- 

Proof. By the assumption and Proposition 5.4, we have 

(:^[J,^",g") ^ (:^o,'^,0) in X D([0,oo),$')2 as n^oo. 

Then, since "R." = ^ ^■^{IZq , S"^ T + and ^ b'^ is continuous by Theorem 3.2 and Proposition 
3.4, the result follows from the continuous mapping theorem with ^'^tc and the addition map. The 
addition map is convergence preserving here by Theorem 4.1 of [32] since the limit of as n — > oo 
is continuous. □ 



6 Diffusion Limits 

We now move on to the diffusion limits. First we define generalized <I>'-valued Wiener process 
and generalized $ -valued Ornstein-Uhlenbeck process as in [3]. These notions will be used to 
characterize our diffusion limits for the age and residual processes. 

Definition 6.1. A continuous -valued Gaussian process W = {Wt)t>o is called a generalized 
$ -valued Wiener process with covariance functional K{s,ip;t,'i/j) = K[{Ws,^){Wt,ip)] if it 
has continuous trajectories and for each s,t>0 and ip,ip G K{s,(p;t,ilj) has the form 

/•sAt 

K{s,(f;t,'ilj) = {Quf,tp)du, 
Jo 

where the operators Qu : $ ^ u> 0, have the properties: 

1. Qu is linear, continuous, symmetric and positive for each u > 0, and 

2. the function u ^ {Quf, V') in D for each ip,ip G 

If Qu does not depend on u>0, then the process is a ^ -valued Wiener process. 

Definition 6.2. A ^' -valued process X = {Xt)t>o is called a (generalized) ^ -valued Ornstein- 
Uhlenbeck process if for each G $ and t >0, 

{Xt,ip) = {Xo,ip) + [ {Xu,Aip)du + {Wt,ip), 
Jo 

where W = {Wt)t>o is a (generalized) ^' -valued Wiener process and A : ^ ^ ^ is a continuous 
operator. 
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6.1 Ages 

Define 

A"" = V^iA"" - A) , A'S = y/^ {A}^ - , = ^ - ^) , = y/^V", pO'" = V^pO'", 

and £^ = -E"5o for n > 1. Then, centering the system equation (14) by the fluid hmit of Theorem 
5.2, for n > 1 we have 

= , - pO' " - ) . 

We now use the following result to approximate V^'^ + by a process that is independent of 
{Aq,E'^) for each n > 1. This is used to prove the required joint convergence in Proposition 6.4. 

Lemma 6.3. Define Z)0'" = (I>°'")t>o G D{[0,oo),D) so that 



"i? / rmAtA{y+n) \ 



1=1 

and c/e/ine i)" = (I)r)t>o G D{[0,oo),D) so that 
/or t > 0, y > 0, where 

-r" = inf {s|nAo(s) > i} , 
= inf \^s\nEs > i} , 

/or i > 0. Zei and 6e t/ie L>([0, cx)), ^>')-?;a/Med processes associated with and D^, 

respectively, for n > 1. Then, if {AqjE"") =^ {Ao,E) in <I>^ x D as n ^ oo, then 

J^pO,n^pn^pO,n^pn^ ^ J^pO^p^pO^p^ !)( [0, oo) , )2 as n ^ CX), (49) 

where T)^ + 1) is a generalized ^'jj^-valued Wiener process with covariance functional given for each 
(p,tp £ o-'iT'd s,t > by 

Kf,o:'p{s,(p;t,ij)= {Auh,(pip) du. (50) 
Jo 

Proof. See appendix. □ 
Proposition 6.4. // iA^,E'^) =^ {Ao,E) in<^_^x D as n ^ CO, then 

+ ^ + in <^'^x D{[0,oo),^'f as n ^ oo, (51) 

where + T) is given in Proposition 6.3 and is independent of E. 
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Proof. By Lemma 6.3 we have the joint convergence 

+ ^ + m X i?([0,oo),$')2 as n ^ oo (52) 

since each of the component processes in the prehmit above are independent of each other. Also 
by Lemma 6.3 we have, 

+ - ^ (0,0,0) in X L>([0,oo),$')2 as n ^ oo. 

(53) 

Combining (53) with (52) gives us our result. □ 
We then have 

Theorem 6.5. // {A}^, S") (Aq, E) in<^>^x D as n ^ oo, then 

A"' =^ A in D{[0, co), $') as n ^ oo, 
where A satisfies the stochastic integral equation 

{At,ip) = {Ao,ip) + EMO) - {V^ + 'D,^) - [ {As,hip)ds+ [ {As,ip')ds. (54) 

Jo Jo 

for all if G //, in addition, E is Brownian motion with diffusion coefficient a, then A is a 

generalized ^_^-valued Ornstein- Uhlenbeck process driven by a generalized ^_^-valued Wiener process 
with covariance functional 

i-sAt 

^i-ivo+v) V') = y (o-^f^o + Auh, Lpip) du. (55) 

Proof. Since = {Aq , —V^'^ —T)^) and ^'^^ is continuous by Theorem 3.2 and Proposition 
3.3, the convergence follows from Proposition 6.4 and the continuous mapping theorem with 
and the addition map. The subtraction map is convergence preserving here by Theorem 4.1 of [32] 
since the limits of P"'" and as n — > oo are continuous. 

If E is Brownian motion with diffusion coefficient a, then i5 is a generalized <l>_^- valued Wiener 
process with covariance functional K£{s,ip;t,ip) = a'^{s A t)ip(0)ip{0). Combining this with (50) 
and the fact that P° + P and E are independent from Proposition 6.4 gives us (55). Thus, ^ is a 
generalized valued Ornstein-Uhlenbeck process. □ 

6.2 Residuals 

Define 

= (^" - ^) , 7t;j = {n^ -Tz), = ^g"", 

for n > 1. Centering the system equation (23) by the fluid limit of Theorem 5.5 we have for all 
n > 1, 
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Proposition 6.6. If {TZ'^.E'^) {TZo,E) in^j^x D as n — > cxo, then 

(7t[},#",a") ^ in <^'tiX D{[0,oo),<l>'f as n^oo, (56) 

where Q is a ^'^-valued Wiener process with covariance functional 

K^{s,^;t,iP) = iE,AEt)Cov{ip{l]),^P{l])), (57) 
where rj is a random variable with cdf F. 

Proof. Notice that for each tp £ (^", ip) can be written as 

Ln, 



By the functional central limit theorem (see Theorem 16.1 of [2]), the first term above converges to 
a Brownian motion with diffusion coefficient y^Var((^(r/)). Thus, our result follows by the random 
time-change theorem (see [11]). □ 

We then have 

Theorem 6.7. // {1Z'^,E'^) (7^o, E) in <^'^ x D as n ^ oo, then 

Tt^^TZ in L»([0, cx)),^>') as n^oo, 
where TZ satisfies the stochastic integral equation 

{TZt, if) = {TZo, ^) + {Qt, + Et{r, [ (R-s, ds, (58) 

Jo 

for all if G $7^. If, in addition, E is Brownian motion with diffusion coefficient a, then TZ is a 
generalized ^^-valued Ornstein-Uhlenbeck process driven by a generalized ^^-valued Wiener process 
with covariance functional 

^s^+g(s, ^\ t, iP) = a\s A t)E[99(7?)]E[V^(7?)] + {E^ A Et) Cov(^(r?), V'(r?)), (59) 

where rj is a random variable with cdf F. 

Proof. Since TZ^ = '^^■r{TZq,£^J^ + G"'), and ^^tc is continuous by Theorem 3.2 and Proposition 
3.4, the convergence follows from Proposition 6.6 and the continuous mapping theorem with 
and the addition map. The subtraction map is convergence preserving here by Theorem 4.1 of [32] 
since the limit of ^ as n ^ co is continuous. 

If E is Brownian motion with diffusion coefficient cr, it is easily checked that ET is a <l> -^-valued 
Wiener process with covariance functional 

K^^is, ip; t, i;) = a\s A t)E[ip{7])]E[^b{rj)]. (60) 

Combining (60) with (57) gives us (59). □ 

Remark 6.8. Notice that in the special case when the arrival process to the nth system is Poisson 
with rate An, E = X- and E is Brownian motion with diffusion coefficient A. Thus, K^j^(s, (p; t, ip) = 
\{s A t)K[{p{rj)ilj{rj)] and Theorem 6.7 gives us a version of Theorem 3 of [8]. 
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7 Markov Process Results 



In this section we prove that under extra conditions on the arrival process and the fluid limit of 
the initial conditions, the limiting age process A of Theorem 6.5 is a time-homogeneous Markov 
process. We then identify the generator of A, which enables us to determine its transient and 
stationary distributions using results from Markov process theory (see [11]). One could also follow 
the same program to analyze the diffusion limit of the residual-service time process in Theorem 
6.7. 

We begin with the following result about the stationary solution to the fluid equation (44) when 
E = X-. This assumption holds, for example, when the arrival process is a renewal process (as is 
the case for that GI /GI /oo queue). 

Proposition 7.1. If E = X-, then A. = XTe is a stationary solution to the fluid equation (44). 
Proof. Plugging A. = AJTg and E = X- into (44) we see that it suffices to verify that 



A / ^{y) dFeiy) = X / ^{y) dF,{y) + At(^(0) - Xt / h{y)^{y) - ^ (y) dF,{y). 

But this follows since 

Xt I h{y)ipiy) - if' (y) dFe{y) = Xt [ h{y)^{y) - ^ {y)F{y) dy 

= Xt[ f{yMy)-F{y)^{y)dy 

= -Xt I {FiyMy))' dy 
= At93(0). 



□ 



The next proposition shows that A has a simpler form when E = X-., E \s Brownian motion 
and ^0 = XJ^e- For each 99 G <1>_4, define the function : <I>_4 — > C by 

F^{^i) = e'^^''^^ for ^Gcl,;^, 

and define the set E{(^X) to be the smallest algebra containing the set {F^p, ip G $^}. We will use 
E{^a) to determine the generator of A. Let C6(<^_4,C) denote the space of bounded continuous 
functions from to C. We then have the following result. 

Lemma 7.2. E{^X) is dense in Cfc($^,C). 

Proof. The result follows by an application of the Stone- Weierstrass theorem for complex-valued 
functions (see Theorem 4.51 of [13]). □ 

Proposition 7.3. If E = X-, E is Brownian motion with diffusion coefficient a, and Aq = XT^, 
then A is a ^ j^-valued Ornstein-Uhlenbeck process driven by a ^_^-valued Wiener process with 
covariance functional 



(73o+©)(s,y';i» = (s At)(o- do + XT^Lptl)). (61) 
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Furthermore, A is a Markov process with generator satisfying 

(GaF^) (^) = ^' - hip) - ^^\0) - ^(^, ] F^i^i), (62) 



for each ip E <^_4, fi G 

Proof. Since E is Brownian motion with diffusion coefficient a, the covariance functional of £ is 
given by 

Kg{s,ip;t,i;) = (s At){a^6o,^i;). (63) 

We now show 

Kpo+^(s,(^;t,V) = A(sAt)(^,(^V^). (64) 

By Proposition 7.1, since = AJTg, A. = \Te solves the fluid equation (44). Therefore, by 
Proposition 6.4, for each ip,ip ^ <^_4, we have 

Kr)0:T,{s,(p;t,tl;) = / (Au,<pi^h) du 







sM 

ip{x)ip{x)h{x) dAu{x) du 







A(sAt) / (p{x)'4j{x)h{x)F{x)dx 
X{sAt) / (p{x)ip{x)f{x)dx. 



(61) then follows by combining (63) and (64) and using the fact that £ and T)^ + !) are independent 
by Proposition 6.4. 

Since Eip{0) is a Brownian motion with infinitesimal variance it^(^(0)^ and {V^ + T), if) is a 
Brownian motion with infinitesimal variance \{J-,tp^), it follows that {At,<p>) is a semimartingale 
and hence by Ito's formula we have for each ip G <I>_4 and t > 0, 

^i{Au^) ^ ^i{A,,^) ^ -^(0) /* e^<-^-^> d4 - i f e*<-^-^> + P„ v?) + i /"* e^<^-^>(A, - h^) ds 

Jo Jo Jo 

-—ip'^{0)[ e*<-^-^> - - / e'^^^'^\j^,ip^)ds. (65) 
Now, plugging in (62) we can write (65) for each ip G and t > 0, 

F^{At) - F^iAo) - [ gAF^U)ds = [ F^{A)d{£s - P° - Vs, ip). (66) 
Jo Jo 

Since ij—P'^—P is a martingale, the stochastic integral on the right-hand side of (66) is a martingale. 
Thus the expression on the left-hand side of (66) is a martingale for each ip G <&^. Since every 
element of £'($^) is a linear combination of elements of {F^p, p G ^a} a^^^d (66) holds for every 
^0 £ ^yi) Lemma 7.2 then implies that A satisfies the martingale problem for Qj[. Applying 
Theorem 4.4.1 of [11] then gives us our result. □ 
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We now wish to calculate the stationary distribution of A. Assuming A has a stationary 
distribution, denote it by vr^ and let denote a random variable with distribution vr^. Recall 
that by the basic adjoint relationship (see Proposition 4.9.2 of [11]) the stationary distribution vr^ 
of A is uniquely determined by the equations, 

a^i^<,(/")vr_4(rf^) = 0, (67) 
for £ Ei<^A)- 

Proposition 7.4. Under the assumptions of Proposition 7.3, Aoo is a ^'^-valued Gaussian random 
variable with mean and covariance functional given by 

E[{A^,ip){Aoo,i^)] = {J'e^iXF + a^F) (^V>, 

for ip,Tp e ^A- 

Proof. By Proposition 7.3, Lemma 7.2 and (67) it suffices to show that for each (p G 

ij^^ F^{fi){fi,^' -hip)7T^{dfi) = (^^^\0) + ^{T,^')^ j^^ F^{fi)7r^{dfi). (68) 

Notice the left hand side of (68) is of the form Kie^^Y, where {X, Y) is bivariate normal with mean 
(0,0) and covariance matrix T, = {T,ij)ij=i^2 with 

Ell = /" (p{xf {XF{x)F{x) +a^F{xf) dx (69) 
Jr+ 

Si2 = ^21= I V{x) (y?'(x) - h{x)ip{x)') {XF{x)Fix) + a'^F{xf) dx (70) 



S22= [ i^ixf {XFix)F{x)+a^F{xf) dx. 



Thus, if we denote the characteristic function of this bivariate distribution by (/)(t) = Ee''^'^^''^'^ = 
e~2* then we can write the left hand side of (68) as 

Eie^^y = 0) = -Si2e"5Sii . (71) 

at 2 

Plugging (69) and (70) into (71) gives us 

i [ _ F^{p){fi,ip' - h<f)Trj^{dn) (72) 
^{x) (^h{x)ip{x) - ip' (x)^ {XF{x)F{x) +a^F{xf) dx j\ F^{fi) ir^idfi). 
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Finally, 

Lfix) (h{x)ip{x) - if'ix)) {XF{x)F{x) +a'^F{xf) dx 



tf^{x)f{x) {XF{x) + a^F{x)) dx - ^{x)^ (x) {\F{x) + a'^F{x)) F{x) dx 



ip^{x)f{x) [\F{x) + a'^F{x)) dx- - f {ip{xf) (AF(x) + (7'^F{x)) F{x) dx 

= ^^'(0) + ^ / ^'(^)f(^) (^^(^) + ^'^(^)) dx + l- [ ^{xf{\ - a^)fix)Fix) dx 

= Y^\0) + ^{J',^'). (73) 

The second equality above follows by integrating the second integral on the left-hand side by parts. 
Combining (72) and (73) gives us (68) and completes the proof. □ 

We can also verify the transient distributions of A using Qj[. Instead of using (67), one may 
use the following generalization (see Proposition 4.9.18 of [11]). Let Ft denote the distribution of 
At for t > 0. Then, (Pt)t>o is uniquely determined by the equations 

/ F^{fi)Ft[dfi]- [ F^{fi)Fo[dfi]= [ [ gAF^{fi)Fs[dfi]ds, (74) 
J^'^ J^', Jo J^'^ 



for F^ G E{^j() and t >0. We then have the following result: 

Proposition 7.5. Under the assumptions of Proposition 7.3, for each t > 0, At is a ^'j^-valued 
Gaussian random variable with mean 

E[(A, = {Ao,F-\^t i^F) > , (75) 

for if G $_4 and covariance functional given by 



K[{At,ip){At,^P)] = \{Te.F~\^t {^^F) {l-F~\^tF)) + / ^(n)V(n) {\F{u) + a^F{n)) F{u)du, 

Jo 

(76) 

for ip,Tp e ^A- 

Proof. We verify the Gaussian distributions defined by (75) and (76) solve the equation (74) for 
each Fip. Plugging the proposed transient distributions into the left-hand side of (74) and using 
the characteristic function of a Gaussian distribution gives us 

Plugging the generator (62) into the right-hand side of (74) gives us 

J^, (i{f^, V>' - V) - Y^'(O) - ^(-F, (^2)^ ds. (78) 
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Next, plugging the proposed transient distributions into (78) and using the characteristic function 
of a bivariate Gaussian distribution as in the proof of Proposition 7.4 gives us 



/ (i (a, F-^t^s [^F -^fyj-\ (^Te, F-\^s [v'F - (1 - Ft.sF) 

J 



(79) 



(^(n) ((^'(2/)F(y) - ^{u)i{^y)) {XF{y) + a^F{y)) du - - ^(•^,'^') 

i{^0,F-ir_,(<pF))-A(^,,F-ir_,(<^^F)(l-F-V_,F))-i/o=<^2(y)(AF(y)+a2F dy 



e 



It now suffices to show that the first factor in the integral of (79) is the derivative of the power of 
the exponential factor. The desired equality (74) will then follow by the fundamental theorem of 
calculus. 

It is easy to see by inspection that the derivative of the first term in the power of the exponential 
in (79) gives us the first term in the first factor of (79). The derivative of the rest of the power of 
the exponential is given by 

-A^ ^[s + y)(if^[s + y)F{s + y)-^{s + y)f{s + y)) - Ek±A^ dy (80) 

+ ^ / ^\s + y)f{s + y)dy- \^\s) {\F{s) + a^F{s)) Hs) 

Using the fundamental theorem of calculus, we can write the last term of (80) as 

- (/.(y) ((/.'F(y) - V9(y)/(y)) (AF(y) + a'^F^y)) dy - "-^^iSi) - - ^ ^^{y)!{y) dy. (81) 

Combining (80) and (81) gives us the rest of the first factor of the integrand of (79) and concludes 
the verification of (74). 

□ 

A Proof of Lemma 6.3 

Proof. We first prove 

+ ^ + P in L>([0,oo),^>') as n^oo. (82) 
By Proposition 4.1 we have 

= n\ ^Jo ^^^'^^'^ ~ du + 2_^ ^{u)'ilj{u)h{u) du 



t rt 



[A^ , ipijh) ds ^ / {As,(p'4'h) ds 
' Jo 
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in D as n — > cxo. The third equality follows from Proposition 2.2 with y > T and (p there set to 
(piph and the convergence follows from Theorem 5.2 and continuity of the integral mapping on D. 
Now by the martingale FCLT, f>^ + T) is a Gaussian martingale with quadratic variation 



«V^ + V»t (99,^) 



(Agh, (pip) ds. 



(83) 



Then, (50) can be proven using (83) and the fact that {V^ + T), 99) has independent increments for 
each if G <I>_4. 

Repeating the argument for P*''" + "D" shows that 



pu,n _^ pn ^ pu _^ P in Z?([0,oo),$ )^ as n^oo. 

It remains to show that the joint convergence (49) holds. To show (49), it is sufficient to show 
that for each ip G <I>^ and t > 0, 



E 



as n — > 00. 



(84) 



This is because since + ^ P° + P and + ^ P° + P individually in L>([0, 00), $') 
as n ^ 00, the sequences (pO'" + P")„>i and (P°'" + P")„>o are both tight, so that ((pO'" + 
P") — (P"'" + P"'))n>o is tight. Then, by Chebyshev's inequality (84) implies convergence of 
finite-dimensional distributions to so that 



(pU,n ^ pn^ _ (pU,n ^ pn^ ^ Q -^^ Z)([0,Oo),$) aS n ^ OO. 

Combining (85) with (82) and using Theorem 11.4.7 of [34] gives us (49). 
We will now show that 



(85) 



E 



Then, the limit 



E 



^O.n \ /^O.n 



as n ^ 00. 



as n — > 00 



can be proven analogously and the previous two limits imply (84). 
Just as in the proof of Proposition 4.1, for t > 0, we have that 



i=l 



where for each n,i > 1 



h[u) du, 



We then have 



(86) 
(87) 



E 



n 



ME 



i=l 



t-T" 



(p{x)d[ l{^,<x} 



ViAx 



h{u) du 
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\nEt\ ^f_fn 




E 

1=1 



/i(n) dti 



t-r" 



Tji Ax 



h{u) du 



< -E ( (^"\ 

\i=\nEt^E];^^ 



\nEt\ 



i=l 



t-fP 



h{u) du 



Focusing on the first term in (88), we have 

2 



1 



n 



[nEtVE^l 
--InEtAEI^I 



-E 



n 



1 

-E 
n 



n 



< -E 
n 



Y: «{vr,^)» 

_i=lnEtAE^] 
■ [uE^yEn pA(t-r,)+ 

/ ip{x)^h{x) dx 

i=[nEiA£t"l ° 

E 

j=[n_EtA£;f 



ip{x)'^h{x) dx 



ip{xfh{x) dxE [\E^ - Et\] 0, 
as n — > oo. Focusing on the second term in (88), we have 



[nEt\ 



i=l 



t-fr 



riiAx 



h{u) du 







-E 



n 



«E 



-E 



n 



LnEtJ 

E « 

i=l 



t-T" 



t-f" 



t-T" 



riihx 



Lp{x) d l{r,,<x} - / Ku) du » 



<p{x)d[ l{n,<x} 



rjiAx 



h{u) du » 
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Now, noticing that 



we have 



n 



< -E 
n 



\nEt\ 



[nEt\ 
i=l 



if{x)'^h{x) dx 



sup_ Irf-f"! 

0<i<lnEt\ 



= inf <^ > - 



n 



-1 / 1 



n 



= inils\E, > - 



sup iTj — Tj 

0<i<lnEt\ 



n 



= sup 

0<i<[n£;t 

< sup 

0<^<Et 



[E-) 



-1 / « 



-1 / « 



n 



E-^ 



< sup I (^'^) \u)-E~\ 

0<u<Et 



U 



0, 



by Theorem 13.7.2 of [34]. Since {supQ<j<^^g^j |r" — 'f"|}„>i is uniformly integrable since it is 
bounded by t, we have proven our result. □ 
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